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SUMMARY 

It is possible to construct migratory dynamo waves in a 
prismatic volume of a conducting fluid, as it exists in the solar 
convective zone. Such waves can explain the principal features of 
the observed solar magnetic activity, particularly the sunspots. 

LIST OF SYMBOLS 

Symbol Definition 

A vector potential 

B magnetic induction 

D electric induction 

e e unit vector in € -direction 

E electric field strength 

f measure for the violence of cyclones 

F Lorentz force 

H magnetic field intensity 

I line integral 

J conduction current density 

k wave number 

K constant 

L separation of periphery from a line C 
p hydrostatic pressure 

q electric charge density 

Q rate of addition of heat per unit mass 
S cross section 


t 


time 


Definition 


Symbol 

T temperature 

U internal energy per unit mass 
v velocity of fluid 

X external force of nonelectromagnetic origin 
r periphery of S 

(6 € £,) Cartesian coordinate system 
€ dielectric constant 

q coefficient of viscosity 

K heat conductivity 

A magnetic diffusivity 

/i magnetic permeability 

p density of fluid 

a electric conductivity 

co angular frequency 

0 expression for angular frequency 

I. INTRODUCTION 

Certain motions of electrically conducting fluids are capable 
of maintaining various magnetic field configurations, such as the 
magnetic fields of the Earth and of the Sun. These mechanisms 
of fluid -magnetic field interaction are called hydromagnetic dynamos. 
This paper demonstrates by means of appropriate wave equations 
that there could exist a plasma fluid-velocity v having a magnetic 
field B which would be maintained by file mutual interaction of magnetic 
field and plasma. The wave equations constitute a hydromagnetic 
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dynamo model. The model is shown to relate to observed sunspot 
formation and activity, and accompanying solar magnetic phenomena. 


II. THE EQUATIONS OF MAGNETOFLUID DYNAMICS 

The behavior of the solar plasma is governed by the equations 
of fluid dynamics and Maxwell's electrodynamic equations. They 
are defined as follows: 

A. Fluid Dynamic Equations 

(1) equation of continuity 

+ div(pv) = 0 

(2) equation of momentum 

+ (v • V) vj = pX + F - Vp 
+ jtjVfdivv) + *lV 2 v (1) 

(3) equation of energy 

P + (v.V)u] = -p div v + K V 2 T + ^ 

+ pQ + q[V 2 v 2 - 2 div(vx curlv) ffcurlv) 2 
- 2 v( VdivT-) - y(divv) 2 ] 


B. Maxwell's Electrodynamic Equations 

curl E = -8B/3t (2) 

curl H = J + 3D/3t (3) 

div B =0 (4) 

div D aq, (5) 

for which constitutive equations are 

B = yH (6) 

D = € E (7) 

J = o(E + v X B). (8) 
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C. Assumptions 

The following simplifications are made so that a practical 
model can be obtained. 

1. The solar plasma is considered as an electrically con- 
ducting single -component fluid. 

2. The magnetic permeability, ju , of photosphere and 
chromosphere is assumed to be that of free space, n Q , and is, 
therefore, a constant (| u = fi 0 = constant). 

3. Viscous forces are neglected (i. e. , the last two terms 
of Equation (1) are set equal to zero), 

4. The displacement current, 0D/3t (in Equation 3), and 
oE (in Equation 8) are neglected, hs well as other electric field 
action on free charges, as in Equation (9). 

D. Deriviation of the Equation of Motion 

Starting with simplified equations, an equation of motion 
is next derived for the convective and magnetic forces which 
dominate action of the hydromagnetic dynamo model. 

The force per unit volume exerted by an electromagnetic 
field on a current, neglecting the action of the electric field on the 
free charges (i.e. , neglecting the term qE),is 

F = J X B. (9) 

In order that the equation of motion for the magnetic field 
may be obtained, E is eliminated between Equations (2) and (8). 
Then Equation (3) is introduced and the identity 

curl curl B * Vdiv B - V*B 
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is applied. When Equations (3) and (5) are used, the equation of 
motion is obtained in the form: 

|f- = curl(vXB) + -i- V 2 B . (10) 

o z o U 

The equation of motion of an electrically conducting fluid 
is obtained through substitution of Equation (9) in the simplified 
equation of motion (1): 

P = p £-f~ + (v . V ) v] = - ^7p + pX + j*XB, (11) 
HI. DYNAMO EFFECTS 

Fluid velocity v and magnetic field intensity H are 
determined by a solution of Equations (4), (6), (10), and (11). The 
question to be answered is whether the conduction current J, pro- 
duced from field B according to Equation (8), can generate the 
field B according to Equation (2) by which J originated. When such 
action does occur it is called a dynamo process. By making the 
following substitutions of terms in Equation (10) 



curl A = B . 


Equation (10) may be expressed in the form: 

= v X curl K + A V 2 A . (12) 

According to a theorem by Cowling, a magnetic field 
symmetric about an axis cannot be maintained by a fluid motion 
symmetric about the same axis. To show this, let C be a line 
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on which B = 0 and about which the lines of induction form closed 
loops which encircle C in the same direction everywhere. Let 
S be a meridian section of C whose periphery r is a closed line 
of induction surrounding C, and let L denote the maximum separation 
of r from C. When Stokes’ theorem and Equations (3) and (8) are 
applied to T , the following result is obtained: 

fa • d? - f curl B • dS = p, JjdS = ^ /(vXB)dS, 

therefore , 

v 

$ B • d?<-^~ f BdS. (13) 

In Equation (13) v max denotes the maximum value of v anywhere 
in the solar volume, and the incremental area dS, lying between 
the lines of induction L and L + dL, cannot in the limit exceed 
d F dL in area. If I(L) denotes the line integral of B over r , it 
follows from Equation (13) that 

/ L I(l)dl. 

X 0 

The line integral 1(b) is a continuous function of the maximum 
separation L, so that I tends to zero as L tends to zero. Let 1^ 
denote the point of the interval (0, L) where I attains its maximum 
value; when Equation (13) is applied to L 0 , the result is: 

l(L 0 )<=^22* L 0 i(l 0 ). 

If I(L 0 ) * 0, the following relationship is obtained: 


X 

V 

max 


L, 
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But L is arbitrary, and if it is assigned a value less than X /v^_ v , 

max 

a contradiction has been established. Therefore, the induction field 
B cannot be maintained in the neighborhood of C. 

For a complete solution of the dynamo problem it would 
be necessary to solve the fluid dynamic and electromagnetic equations 
simultaneously. Because the difficulties of a complete solution 
are too great, the more restricted problem is usually investigated 
to establish whether a motion exists which can maintain a magnetic 
field. 

A magnetic field which is everywhere perpendicular to 
meridional planes is called "toroidal, " and a magnetic field whose 
field lines are everywhere in meridional planes is called "poloidal. " 

A rectangular Cartesian coordinate system (6, €, £) will be 
applied to the fluid. The vector potential A is considered to extend 
in the € -direction. Writing A for the € -component of A, we can 
express Equation (12) in the form: 

|4 = fB + A V Z A. (14) 

o t 

It will be assumed that a moving fluid rotates about the 
j^-axis. Dynamically, such a fluid motion is related to the cyclones 
and anticyclones observed in the Earth's atmosphere. Such a 
cyclonic motion results when the Coriolis force component adds 
to the primary driving force of radial convection. Similarly, 
deflection along the length of a vertically rising and diverging 
column of solar plasma would tend to produce a cyclonic configuration. 
The factor f in Equation (14) is a measure of the violence of such 
cyclonic motions. 
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If fluid motion is present: 


v = e € v; 


then, representing the nonuniform rotation by Equation (10), the 
toroidal field B is generated as a function of the magnetic field 
potential A: 


|f- = V X [ v X ( V XA) ] + A V*B. 
d t 

Since A and v are independent of € and have only € - 

d § 

components, -r*— is in the C -direction, and the last equation 
o ^ 

reduces to: 


3B / 3v 3 A s v 3 A \ 

at = \ a; a 6 * 36 3 ^ / 


+ XV 2 B. 


(15) 


In order to solve the dynamo Equations (14) and (15), we assume 
uniform shearing. If v varies linearly with £ : 


= K = constant, and 


= 0 . 


Then Equation (15) reduces to the form: 


If' K 17 (16) 

If the cyclones are distributed uniformly throughout the 
space, K is a constant, and the dynamo equations reduce to two 
simultaneous linear equations with constant coefficients. 

The following trial solutions are made regarding the 
solutions of Equations (14) and (16): 

A=A 0 e l <" t + W) 

B s B 0 e i < wt + k6 > . (17) 
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After substitution of Equations (17) into Equations (14) and (16), 
the following results are obtained: 

A 0 = (i w + Xk?)*B 0 f=O 

and 

\ 

-A 0 ikK + B 0 (i w + Xk 2 ) = 0. 

Setting the determinant of the coefficients equal to zero yields the 
result: 

(iu + Xk 2 ) 2 - ikKf = 0, 
which can be expressed in the form: 

1 

iw + Xk 2 = ± (ikKf) 2 = 0(l±i); 


therefore. 




•(^)* • 


ico = ( S2 - X k* ) ± i SI , 


and 


A 0 = -B 0 n 


(&)■ 


We now use t’ e last two equations with the trial solutions (Equations 17) 

of the dynamo equations (Equations 14 and 15) to obtain equations 

describing the migratory dynamo waves: 

„ _ (SX-Xk 2 )t i(k6 ±€lt) 

o — o e e 

o 


A = 


JL 

kK 


(l±i)e (n “ Ak)t e i(k6 ±0t) 


The solar dynamo, if it exists, would be located in the convective 
zone in a shell sufficiently thin ( « 1C? km) that its curvature may be 
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neglected. Application of the dynamo equations to such a flattened 
space yields just the migratory dynamo. 

Dynamo waves migrate from the polar to the equatorial 
regions in the convective zone just below the surface of the Sun. 
Because the liquid -dense photosphere has a relatively high sonic 
and magnetic conductivity ir. comparison with the adjacent gaseous 
chromosphere, the dynamo waves are constrained from diffusing 
into the chromosphere. Thus only in regions of most intense 
cyclonic motions or of nonuniform rotation will significant wave 
strands be carried beyond the normal radial limit of the photosphere 
It is in such regions that sunspot activity would be observed. 

The Sun's magnetic field is composed of traveling magnetic 
waves. The traveling wave is assumed to consist of two components 
toroidal and poloidal, as indicated diagramatically in Figure 1. 

The sense of the latitudinal toroidal bands alternates from one band 
to the next. The poloidal field components are in meridional planes 
rotated 90 degrees from the planes of the toroidal components. 

In each hemisphere, there are two or three toroidal bands 
at any one time. About 22 years are required for migration of a 
band from pole to equator, so that new bands approach the equator 
at 11 -year intervals. The poloidal component of the traveling wave 
predominates from pole to middle latitudes, the toroidal component 
dominates below the middle latitudes, and both components vanish 
at the equator. 

An occasional strand of plasma, as a result of particularly 
intense local wave amplification, moves to or above the surface 
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of the photosphere, but the main field configuration remains 
essentially unaffected. Sunspot activity peaks at 11 -year intervals, 
as new toroidal bands in each hemisphere move into the region in 
the lower latitudes of most intense field amplification. Strands 
thrown out below the middle latitudes are predominantly a result 
of action of local toroidal field intensification, while strands 
occurring nearer the pole result from local poloidal field buildup. 

Below the middle latitudes the spots occur in pairs. These 
paired spots correspond to the exit and entrance locations of in- 
verted U-shaped distortions of the toroidal flux tubes which have 
been pushed above the normal radial limit of the photosphere. 

The east-west orientation of the sunspots in the lower 
latitudes results from the east-west direction of the dominant 
toroidal component of the flux tube. The migration of a spot toward 
the equator from its region of formation is a result of the migration 
of the associated toroidal field. The reversal of polarity of the paired 
spots corresponds to the alternation of field each half cycle (that 
is, with each successive toroidal band). 

IV. RESULTS AND CONCLUSIONS 

The explanation of the Sun' s magnetic field given in this 
paper is derived from studies by E. N. Parker, and may be 
summarized in the following way. 

The Sun's magnetic field is generated by, and moves in 
unison with, the electrically conducting solar plasma. Hence, 
fields of major intensity are confined to the convective zone. 
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Field and moving plasma continuously interact to modify their 
interdependent configurations. The field is azimuthal at low 
altitudes, and progressively changes toward a meridional plane 
near the poles. 

The field is maintained as a result of nonuniform axial 
rotation of successive layers of the Sun's fluid mass and 
vortex effects in the convecting plasma. The nonuniform rotation 
amplifies the azimuthal component of the force field. Vortex 
effects in the convecting plasma produce field components in 
poloidal planes. The net effect is to twist the azimuthal field 
loop by a process of reenforcing the loop on one side and reducing 
field density on the other. Simultaneously, the dynamo effects 
intensify the field. Similarly, the magnetic force-field differential 
causes surface loops to move toward the equator and subsurface 
loops to move toward the pole while dynamo effects act to intensify 
the fields. 
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FIGURE 1. Diagram of Solar Toroidal and Poloidal Magnetic Fields 
in the Outer Half of the Convective Zone (from Parker, Ref. 1) 


